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Abstract 

In this paper, we are interested in the analytical study of a nonlinear 
Stochastic Partial Differential Equation (SPDE) arising as a model of 
phytoplankton aggregation. This SPDE consists in a diffusion equation 
with a chemotaxis term responsible of self-attraction of phytoplankton 
cells and a multiplicative branching noise. Existence of mild solutions is 
established through weak and tightness arguments. 
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Introduction 

In this paper, we are interested in the following nonlinear stochastic partial 
differential equation: 

g d ' 2 d _ • 

—w(t, x) = D—u(t, x) - — (■ u(t , x) [G * u°(t ,.)] (a;)) + y/\u+(t,x)W{t, x), 

(1 ) 

in [0,T] x C, where Q = ]0,T[ is a bounded domain with boundary dfl in 
R, a; is a one dimensional coordinate, t is time. The motivation in studying 
this SPDE is that equation m arises as a model of phytoplankton aggregation 
(EimillMEI). In fact, the authors in ([9j, mm .fT2l) have investigated an 
individual-based model (IBM) of a population of phytoplankton that takes into 
account small scales biological mechanisms for phytoplankton cells. These pro¬ 
cesses are: random dispersal of phytoplankton cells due to turbulence, spatial 
interactions between phytoplankton cells caused by chemical signals and ran¬ 
dom division and death of phytoplankton cells. The aim of such modelling was 
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to capture the main features of the individual dynamics at small scales that 
are responsible at a larger scale for the formation of aggregating patterns. An 
Eulerian version of the IBM has been derived in (isi. mi) and the SPDE <[T|) is 
obtained as a continuum limit of this IBM. rt(f, x ) represents the spatio-temporal 
distribution density of phytoplankton on the vertical water column. We denote 
by u+(t,x) the positive part of u(t,x ), that is u+(t,x) = max(it(i, x), 0). The 
diffusion term in m takes into account the random spread of the phytoplank¬ 
ton cells with the coefficient of diffusion D; while the transport term, which is 
a chemotaxis term, describes the interaction mechanisms between the phyto¬ 
plankton cells via the velocity G *u°. The latter has the form of convolution as 
in 22], i.e., 

[G * u°(t,.)] (x) = f G (x — x') u°(t, x')dx r , 

Jr 

where G is the attractive force given by: 


and 


G (x) [( 1^1 fl) ( 1^1 r o)] 1 [— n, — r- 0 ]U[ro,ri] i x ) • 


f u(x) 0 < x < L 

\ 0 x < 0 or x > L. 


The biological explanation of the interactions between phytoplankton cells is 
based on the following works ( [H] . [ 13 ] . [ 15 ] . pH ] . [22 ] . [ 21 ] . [ 25 ] . [ 25 ] ) which report 
that species of phytoplankton such as dinoflagellates and motile algae have 
chemosensory abilities i.e., they can sense the chemical field generated by the 
presence of other dinoflagellates and motile algae that are present at a certain 
distance. More precisely, dinoflagellates and motile algae leak organic matter 
such as sugar and amino acids, forming regions around them having concentra¬ 
tions higher than average |3j. Experiments studies on the chemosensory abilities 
in dinoflagellates and motile algae (ee mi) show that the released products 
attract other dinoflagellates and motile algae that are present in a suitable neigh¬ 
borhood. It has also been observed that high concentrations of these products 
inhibit the chemosensory behavior in dinoflagellates and motile algae EE 
So, if a phytoplankton cell is located at a position x, the extracellular products 
released by this cell form a concentration around x, which is highest in the 
closest vicinity of x (for instance, on a radius of length ro ( ?'o small)) and then 
decreases progressively. Via their chemosensory abilities, all dinoflagellates at 
positions x' such that ro < \x — x'\ < n ( r± > ro and ro too small relative to 
ri) detect the differences of concentration in water and hence are attracted to 
the cell in x. Beyond n, cells cannot perceive the difference in concentration 
because they are sensory limited (IE EE EE- Hence, they are not attracted. 
(G * u°) (x) describes the velocity induced at the site x by the net effects of 
all phytoplankton cells at various sites x'. The kernel G{x — x') associates a 
strength of interaction per unit density as a function of the distance x — x' be¬ 
tween any two phytoplankton cells in sites x and x'. G behaves as a gradient, 
that is, a cell is attracted to the region of high density. 

W(t,x) is a white noise in space and time defined on some probability space 
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(AjJ 7 , P) j27] and A is the branching rate. The term y / Au + (t, x)W(t , x) de¬ 
scribes the stochastic fluctuations of the number of phytoplankton cells as a 
result of the random birth and death events. 

Here, R represents the vertical axis oriented downward from the surface to the 
seabed. The point 0 is at the surface of water and L is the limit of the ”eu- 
photic zone” ie. the upper layers in the vertical water column. The restriction of 
model © to this domain is due to the fact that phytoplankton cells can survive 
and multiply only in the ’’euphotic zone”. Therefore, von-Neumann boundary 
conditions are imposed at the surface 0 and at L: 

d 

—u(t,x) = 0, on [0,T] x dfl. (2) 

and the initial density is 

u{0,x) = uq(x) > 0, in Cl. (3) 


To our knowledge, SPDE © is new and unknown in both biological and 
mathematical literature. From the mathematical point of view, it is a complica¬ 
tion. The major difficulties come from the two nonlinear terms: the chemotaxis 
term with a convolution and the stochastic branching term with the nonlinear 
multiplicative noise. We point out that previous works ( 0 , 1 ©]) have concerned 
the study of the deterministic part of the SPDE © (i.e. equation © without 
the noise term) with the goal of analyzing the influence of the spatial inter¬ 
actions between phytoplankton cells on the aggregation process. Here in this 
paper, we are interested by the whole equation © which includes the effect 
of the stochastic branching term. The volatility \J Ait+(t, x) of the space-time 
white noise is non-Lipschitzian; hence, we are not able to expect the existence 
of strong or mild solutions for ©■ However, we might expect the existence of a 
mild solution to © in the weak sense; that is we shall construct a probability 
basis (A, J 7 , (^Otefo.T]:-P )j on which there exists a Gaussian space-time white 


noise W(t,x) and a mild solution u(t,x ) to ©. This will be done by using 
tightness arguments which are useful in obtaining weak convergence. The rest 
of the paper is organized as follows. In the next section, we present the abstract 
formulation of the problem and in section 3, we present our results on existence 
of mild solutions at the weak sense. Several steps are arranged for proving the 
main theorem. Some estimates on the semi-group generated by the operator 
d 2 

—^- are listed in the Appendix. 
ox z 


1 Abstract formulation 

We reformulate problem © with conditions © and m as a stochastic version 
of an abstract Cauchy problem which can be treated by using the theory of 
semigroups of operators: 

/ du(t) = (Au(t) - B[u(t)g G (u(t))])dt + C(u(t))dW(t) 

\ u(0) =u 0 W 
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The operator A : T>(A) Cl : = L 2 (fi) — > X is defined by 


Aw 

V(A) 


= D- 


dx 2 ’ 

= lweH 2 (n ): 


'| an 


= 0 


}■ 


and the operator B : V{B ) C X —> X by 


( 5 ) 


Bw = —to, (6) 

dx 

V{B) = 

H 1 ( fi) and H 2 (tt) denote usual Sobolev functions spaces. We will denote by (,} 
and ||.|| , respectively, the scalar product and the norm in X. The operator A 
commutes with B and they are related by the following formula 

(Bu, DBu) = - (u, Au ), Vm G V(A). (7) 

We endow D(B) with the graph norm \x\ v ^ = ||Bx|| + ||x|| for x G D(B). 
The operator go is defined as follows: 


9G (y>) (x) = [G * ip\ (x) 


G(x — y)ip {y) dy. 


( 8 ) 


By straightforward consequence of standard calculations, we can establish that 
go ■ T>(B) —» T>(B), continuously so there exists a constant 5, so that 

I 9 g ^ (9) 

Note also that G*u° is uniformly bounded. Hence, gdd is uniformly bounded. 
As a result of Holder’s inequality, we get 

vTIGUIHI, Vue D(B). (10) 

On the other hand, we have 

l%WL < \u \ v(B) , Vu e D(B). (11) 

C is the non linear operator, 

C:X ^L hs (V(B),X) 

where C(u) is the linear operator of multiplication by the function y 7 A u +, that 
is: 

Vu G X, w e B(B), (C(u) w)(x) = y 7 A u+(x)w(x), x e 
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Lhs(P{B),X) denotes the space of Hilbert-Schmidt operators from V(B) to 
X equipped with the norm 


llTlI^^HTefell 2 ) 1 / 2 < oo, 

k =1 

where {e*,} is a complete orthonormal basis in X. 

Suppose that the compact embedding: 

J : T){B) > X is Hilbert-Schmidt, (12) 

then (W(f)) te [o,T] can be considered as a cylindrical Wiener process on V{B) 
with values in X , defined on the probability space (A, J-,P ) with a filtration 
t)te[o,T] ■ 

Proposition 1 The operator A defined by (0 is the generator of an analytic 
semigroup of contractions in X, (T(t))t> o, compact for t > 0. The restrictions 
T{t)/D(B) send P(B) into itself and are uniformly bounded in T){B) (that is, 
there exists C\ > 0, such that, \T(t)/x>(B)\ b(b) — @ 1 , for t > 0). 

Proof. The proof is similar to that one in [2]. ■ 

We propose to solve, on time interval [0,T], equation (HJ) in integrated form 
by using the stochastic generalization of the classical variation of constants 
formula 

u(t) = T(t)uo— f T(t — s)B[u(s)gG (u(s))]ds + f T(t — s)C (u(s)) dW(s). 
Jo Jo 

(13) 

We remind that a predictable X-valued stochastic process (u(t)) te i 0 ,t] defined 
on the probability space (A, T, (J c ’t)te[o,r]i-f > ) is called a mild solution of the 
differential equation if u(t) is a solution of (fl3l) . On another hand, the 
stochastic process (w(i))te[o,T] is called a weakened solution of ® if u(t) is a 
solution of 


l Z Z 

u(t) = uo + A(J u(s)ds) — J B [u(s)gG (u(s))] ds + J C(u(s))dW(s). (14) 

oo o 

It is well known that solutions of (fl3l) and m are equivalent ABLE)- Note 
that the stochastic integral / T(t — s)C (u(s)) dW(s) is well defined since 

t ^ 0 

E{ / \\T(t — s)C (u(s ))||2 ds) < oo ,Vt S [0, T] . This is due to the fact that 

J o 

C (u(s)) is Hilbert-Schmidt operator and T(t — s ) is linear bounded operator 
then, basing on the theory of Hilbert-Schmidt operators (e.g. |TQ, Chapter l), 
T(t— s)C (u(s)) is Hilbert-Schmidt operator too. 
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2 Existence of solutions 


This section is concerned with the existence of mild solutions for ©■ For this 
purpose, we need first to establish several propositions and lemmas. We start 
by giving some useful estimates. 

Lemma 2 1) There exists a constant M, such that, for all u, v G 25(1?), we 
have 

II B[ug G {u)} -B[vg G (v)]\\ < Mmaxflul^ , \v\ v(B) )\u - v\ v(B} . 

2) There exists a positive constant Q, such that, for all u £ V ( B ), it holds that 

\\B [ug G {u)}\\ <Q\u\ v(b) \\u\\. 

3) There exists a positive constant C , such that, for all u G X, it holds that 

\\BT(t)u\\ < |M|, Vf > 0. (15) 

Proof. For the claim, the proof is the same as the one given in (Lemma 2.1, 

OP)- ■ 

By using Lemma 2, we can prove 

Proposition 3 Consider the following abstract problem 

f du(t ) = (Au(t) - B [u(t)gc (u(t))]) dt + a(u(t))dW(t) 

\ u( 0) = u 0 . ^ 

Suppose that the operator a satisfies 

a:X^L HS (V (B),X) 

and there exists a constant I\ such that for all x, y € A', 

l|aO) - a(y)|| 2 < K\\x-y\\. 

Then, for each uo € D(B), problem \16\) has a unique mild solution u(t), t G 

[0,T] such that sup E(\u(t)\ D , B A < oo. 
te[o,T] 

Proof. We use the method of successive approximations. Let uo G D(B) and 
define the sequence (u n ) n > i by 

t t 

u n +i(t) = T(t)uo(t) - J T(t-s)B[u n (s)g G {un(s))\ds + J T(t - s)a(u n (s))dW(s). 
o o 

Let us denote by Y the Banach space of all nonanticipating D(13)-valued con¬ 
tinuous stochastic processes {L/(i)} tg [ 0 T j endowed with the norm 

f 1 1 > 2 

|||t/|||= sup S|t/(t)|^ (B) < OO 

lo <t<T ) 
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Assume that (u n (t)) n > 1 is bounded in D(B), that is, 3 R > 0 such that \(u n (t)\ D ^ < 
R, Vt € [0,T],Vn > 1. Hence ( u n (t),t € [0,T])„>i is also bounded in Y, since 

r 2 1 1/2 

|||u„||| = < sup E \u n {t)\ D(B) > < R. 

Lo <t<T ) 

Let 

hn+i(t) = E(\u n+1 (t) - u n (t) lo(B)), n > 0 

We have 


hn+i(t) < 2 E 


T(t - s ) [a(u n (s)) - a(u„_i(s))] dW{s) 


D(B) 


+2 E 


T(t - s) ( B[u n (s)g G (u n (s ))] - B[u n -i(s)g G (u n -i{s))]) ds 


< 2 E(J || T(t - s ) [a(u n (s)) - a(u„_i(s))]|| 2 ds) 
o 


+2 E 


T(t - s) ( B[u n (s)g G {u n (s ))] - B[u n ^i(s)g G (u n -i(s))]) ds 


ii 


2 

D(B) 


2 

D(B) 

✓ 


Since 3 Ci > 0 such that T(t) /v(B) d ^ < Ci, for t > 0, hence 


I < 2 CfE( J ||a(/u n (s)) - a(u n -i(s ))\\2 ds) 

o 

t 

<2 CfK 2 E(J |M S )-u„_ 1 ( S )|| 2 ( fe) 

0 

< 2Cf K 2 T E( sup |u„(s) - u n -i(s)\ 2 D(B) ). 

0 <s<T 
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To calculate //, we have 


T(t - s ) (B[u n (s)g G (u n (s))] - B[u n - 1 (s)g G (u n -i(s))]) ds 


D(B) 


T(t - s) (B[u n (s)g G (u n (s))] - B[u n -i(s)g G (u n -i(s))]) ds 

t 

B J T(t — s) (B[u n (s)g G (u n (s))] - B[u n -i(s)g G (u n -i(s))]) ds 
o 

t 

< J Mmax(|u„(s)| D(B) , |u„_i(s)| £)(B) ) |it ra (s) — u„_i(s)| £)(B) ds 


C 


\/t - 


\\B[u n (s)g G (u n {s))} - B[u n - 1 (s)g G (u n -i(s))]\\ ds 


<MR(T + 2 VTC) sup \u n (t) — u n -i(t)\ 

o <t<T 


D(B) ■ 


Therefore 


II < 2 M 2 R 2 {T + 2 VTC) 2 E{ sup | u n (t) - u„_i(i)|Lm)- 

0 <t<T 


Hence, we obtain 

h n+1 {t) < 2[M 2 R 2 ( T + 2\/TC) 2 + C 2 K 2 T]E{ sup \u n (t) - u n - X {i)\ 2 ) 

0 <t<T 


By choosing T > 0 small enough so that 

2 [M 2 R 2 ( T + 2 VTC) 2 + CfK 2 T ] < i, (17) 

we obtain 

1 2 

^n+i(i) < sup |ui (t) - u 0 {t)\ D(B) ). 

u n +i(t) — u n (t) is the general term of an absolutely convergent series in the 
Banach space Y. Hence, {u n (.)} is a convergent sequence in Y. The limit in Y, 
Uoo(t) = lim u n {t) is a solution of (fTfil) for fixed t. To complete the proof of 

n —>-oo 

the proposition, we have to show that the sequence u n {t) remains in Y, i.e. 
E( sup \u n+ 1 (t)\ 2 D(B) ) < oo. 

0 <t<T 


E(\u n +i(t)\ 


< 2 C 2 \u 0 {t)\ 2 D(B) + 2 M 2 R 2 {T + 2 CVT ) 2 

E( sup \u n {t)\ 2 D(B) ) + 2ClKlT E{ sup \u n (s)\ 2 D(B) ) 

0 <t<T ' ' 0<s<T v 7 


< 2 ClR 2 + 2 R 2 


M 2 R 2 (T + 2 cVT ) 2 + C 2 K 2 T , Vt G [0, T]. 














Using T > 0 such that 1171) holds, we have 


sup E(\u n+1 (t)\ 2 D(B) ) < 2 ClR 2 + \r 2 < oo . 

0 <t<T * 

To prove the uniqueness of the solution (up to equivalence), consider u and v 
two mild solutions of (fldl) on [0,T]. We might show that 

sup E\u(t)-v(t)\ 2 D{B) =0. 


We have 


E\u(t) -v{t)\ 2 D(B) <2E 


T(t - s)(B[u(s)g G {u{s))} - B[v(s)g G (v{s))])ds 


D(B) 


+2 E 


(i) 


T(t — s)[a(u(s)) — a(u(s))]dW(s) 


D(B) 


We obtain 


( 2 ) 


(2) < 2C(K Z TE( sup |«(a)-«(a)|D (B) ) 

0 <s<T v ' 


To calculate (1), we have 


T(t - s){B[u{s)g G {u{s))} - B[v(s)g G (v(s))])ds 


<MR(T + 2CVT) sup j'u(f) — v{t)\ D , B \ 

o <t<T y 1 


D(B) 


Then 


(1)<2M 2 1 2 R 2 (T + 2CVt) 2 E{ sup \u(t) - v(t)\ 2 D(B) ) 
Finally, we obtain 


0<£<T 


E\u(t)-v(t)f D{B) < \2M 2 R 2 {T+ 2CVT) 2 + 2CiK 2 T 

E( sup \u{t) - v(t)\ 2 D{B) ) Vte[0,T] 


0 <t<T 


Since T is chosen such that 
have 


2 M 2 R 2 (T + 2 CVT) 2 + 2C 2 K 2 T 


< 


sup E(\u(t) 
0 <t<T 




hence we 
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which leads to 

u(t) = v(t), Vi £ [0,T], P — a.s. 

■ 

We return to problem f4]). We propose as in [5] to construct a sequence of 
Lipschitz continuous mappings (a„(«))„gN to approximate the non-Lipschitzian 
function^A m+ on R. For n £ N, define a n (.) by 

{ 0 u < 0 

0 <u<- ( 18 ) 

y/Xu u > — 

n 

It is clear that the sequence (a n (u)) n converges to V A u + uniformly on u £ R 
as n —> +oo. 

Lemma 4 Suppose that the sequence (a n (u)) n defined by \18\) satisfies 

a n :X^L HS (V(B),X) (19) 

so that 

\/u £ X, w £ 'D(B), a n (u) w = a n (u)w. (20) 

Then, there exists a constant K' such that for all u, v £ X, 

II a n (u) - a n (v )|| 2 < K' ||u - u|| . (21) 

Proof. Let u,v £ X and {e^} a complete orthonormal basis in D{B). Hence 

2 + °° 2 
||a n (w) - a n (u)|| 2 = E \\(an(u) - a n (v))e k \\ 

fc=l 

+ 00 ^ 99 

= Z / \a n {u(x)) - a n (v(x))\ \e k (x)\ dx 
k= 1 0 

As a n (u ) is lipschitz in u on R, then there exists I\ > 0 such that 

o +oo L 

II a„(u) - a„(v )|| 2 < K 2 E f \u(x)) - v(x)\ 2 \e k (x)\ dx 

k= 1 0 
+00 

<K 2 E \\{u-v)e k \\ 2 
k= 1 

+00 

<K 2 \\u~v\\ 2 j:\e k \ 2 D{B) 

k—1 

Using the hypothesis m, we have 

+ OO +OO 

l efe l d(b) = II Jefc H 2 = II J ll2 < 00 » 

k= 1 k=l 
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which leads to the fact that there exists M\ > 0 such that 


||a n (u) - a n (v)\\l < M 1 K 2 ||u - v\\ 2 


Let us now consider the following approximating SPDE’s: 


u n {t) = Io(t) - I n {t) + J n (t), neN (22) 

with 

Io(t) = T(t)u 0 (t) 

t 

In(t) = j T(t- s)B[u n (s)g G (u n (s))\ds 

0 

t 

J n (t) = J T(t-s)a n (u n {s))dW (s) 

0 

where the sequence («„(.))„ is defined by (fT51) and satisfies (fTUl) . 

We aim to prove weak convergence of the sequence {u n {t), t £ [0, T]} raS N by 
proving the tightness of {I n (t),t £ [0,T]} n£ N and {J n {t),t £ [0,T]} raS N since 
Io(t) is deterministic. 

Proposition 5 Let the initial condition uq £ D{B ) be a continuous determinis¬ 
tic mapping on fl.The sequences {I n (t),t £ [0,T]}„ e N and {J n (t),t £ [0,T]} ng N 
are tight in C([0, T], D(B)). 

Proof. By Lemma 4, it follows that for each n £ N, the sequence (a n (.)) n 
satisfies the conditions of Proposition [3] Hence, there exists a unique mild 
solution u n (t),t £ [0,T] for problem (B^l) such that sup E{\u n {t)\ 2 D(B ^) < oo. 

te[o,T] 

For given 0 < t' < t < T and rj £ [0, i], it follows from Lemma[T] (in Appendix 
A) that 

t 

E\J n {i) - J n (t')\ 2 D(B) < 2E{ J \\(T{t — s) ~ T{t' — s))a n (u n (s ))||2 ds) 

0 

t' 

+2 E(J \\(T{t' - s)a n {u n (s))\\lds) 

<2C'K 2 E{ sup \u n (s)\ 2 D , B) )\t-t'\ ri 

0 <s<T v 

t' 

+2 E( J \B(T(t' - s)\ 2 d{b) \u n (s)g G (un(s))\ 2 DiB) ds 


ii 


On another hand, by using Lemma [5] (in Appendix A), we have for given 
0 < f < t < T and rj € [0, 


E | I n (t) - I n (t')\ 2 D[B) < 2 E( J | (T(t - s) - T(t' - s))B[u n {s)g G {u n {s))]\ 2 D(B) ds 

0 

t' 

+ 2 E{ J \B(T(t' - s)[u„(s)sg(u„(s))]|^ ( b) 

t 

t 

< 2 E( J \B(T(t - s) - T{H - s))\ 2 d(b) \u n {s)g G (u n {s))\ 2 D(B) ds 

0 

t' 

+2 E{ f \B(T(tf - s)|^ (B) \u n (s)g G (u n (s))\ 2 D{B) ds 


using ©, we obtain 


E\I n (t)-I n (t')\^ B} < 2C"5 2 E( sup \u n {sT D{B) )\t-tr 

v ' 0 <S<T v ' 

+2 C'iS 2 E( sup |ltn(s)|^ (B) ) \t- t'l 71 

0 <S<T v ' 

<c 2 |t-t , r. 


It follows that 
T T 

E 



0 0 


I Jn(t) Jn{t')\D(B) 


T T 


dtdt' < J J Ci\t — t , | 7) 27 dtdt' 

0 0 

< 00 for all 0 < 7 < 


For R > 0, we define A R C A 

T T 

A R = {log A / 

'0 0 



|J„(t,a;) J n {t', w)|£)(b) 


|t — t‘ 


/|7 


dtdt' < R 


If we let in the Garsia-Rodemich-Rumsey Lemma (Lemma 9, in Appendix A): 
/ = </„(., w), x = t, y = f', p(ir - t/) = |ar - y | 7 *e p(x) = |x| 7 and ip(x) = x 2 , 
then 

Ml 

\\Jn{t,uj) - J n (t',Uj\\ DiB) < 8 J 1 r 1 {^;)dp{u). 

0 
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Since ip 1 (^5) = supv = 2 A-A ; thus 


s/R 


2 <- 4 R U 


V 2 < 




\Jn{t,u) - J n {t',u\ D{B) <8 J 2-^7 u 1 1 du 


< 16 — f — i?5 It-t'l 7-1 

7-1 


Finally 


It follows that for any 0 < S < 7 

| — J n [t ? ^|x?(_B) 


Sup 


I *-tf 


< sup Cilt-t'l 7-1 5 1? 1 / 2 


Hence, for w £ Ar, we can define 


ll^ra(-> w )llc < ([0,T],U(B)) — S7l P j l^(^’ W )I_D(B) 

|J n (t,w) — j n {t ^)\r>/r\ 

sup 


< C'R 1 / 2 




We have then for every P > 0 

P(A R ) =p{weA/ 
< P 


T T ''Ut)-J n {t')\ D(B) \ 2 

\t-tV ) 


dtdt' < R 


0 0 


{ll^"(-! W )llc i 5 ([ 0 ,T],D(B)) — ^ ^-} 


Let us now show that {J n (t),t £ [0,T]} ne N is tight on C([0, T], D(B)). For 
R > 0, we define 

Br = {/ € C([0,T],D{B))/\\f\\ c7ao nD{B}) < p) . 

Br is a compact set of C([0, T], D(B)) (see the Proof B1 in Appendix B) and 
we have the following: 


SUpP[Jra £ Br] — SUpP[|| Jn(-, ^ 110 5 ([0 T] ,D(B)) > B] 
n n u ’ J ’ v 


< 


< 


sup P 


supP 



J\n(t) Jn(d']\ D(B) 


0 0 


1 1 — t‘ 


/IT 


dtdt' > 


R 2 
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Using Markov inequality, we obtain 


— c " 2 

sup P[J n G B r \ < — sup E 

n -tb n 



^ ) 



(23) 


Since E 



/ \Jn(t) Jn{t')\u(B) \ 

^ \t^W ) 



< +oo, (E51) leads to 



supP[J ra G Br] < 

n -Lb 


< e 


for R large enough. This proves the tightness of { J„(f), t G [0, T]} ng N. A similar 
argument implies that {I n (t),t G [0,T]} ng N is also tight on the same space. ■ 
Now, we are able to state our main result 


Theorem 6 Let the initial condition uo G D(B) be a continuous deterministic 
mapping on f2. Then, there exists a probability basis (A, T, (.Ft)t g [o,T]> P) on 
which there is a cylindrical Wiener process (W(t)) te ^t] and a mild solution 
i u (t))te[o,T] of OW inC([0,T],D(B)). 

Proof. The tightness of {I n (t),t G [0,T]}„ 6 n and {J n (t),t G [0,T]} ngN implies 
that {u n (t),t G [0,T]}„ g N converges weakly on C([Q,T],D(B)), this means that 
there exists a subsequence {Nk}keN C N and a U([0, T\, D(f3) (-valued random 
variable v such that 


UjV fc —> v on C([ 0 ,T],I 1 (B)), as k —> +oo. 

By Skorohod Representation Theorem, there exists a probability space (A, T , P) 
with a filtration (Pt)te[o,T] and C([0, T], 15(1?) (-valued random variables (uAr)./veN 
and u such that as N —> +oo 


and 


un —} u P-a.s on C([0,T],D(B)) 


^ in law . in law 

un = u n and u = v. 


The mild solution un {-) of (B^l) is also a weakened solution of (B^l) (see for 
instance 0 ), that is un(.) satisfies the following integral equation 

t t t 

u N (t) =u 0 + A(J u N (s)ds)~ J B [u N (s)g G {u N (s))]ds + J a N (u N (s))dW(s). 
oo o 

Hence, for each N, 

t 

M N (t) = J a N (u N (s))dW(s) 
o 

t t 

= u N (t) -Uo- A(j u N (s)ds) + J B[u N (s)g G ( u N (s))]ds 
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is a square integrable martingale with respect to = a {uat(s), s < t} and has 
the following quadratic variation 


(• M N (t),M N (t)) = / a N (u N (s))a* N (u N (s))ds. 


(24) 


Since £{un) = £{un) {£ denotes the probability distribution), 


M N (t ) = u N (t) -u 0 -A( u N (s)ds) + / B [u N (s)g G (u;v(s))] ds 


has the same distribution as Mjv(t) and hence 


E 


MN (t) = E\M N (t)\ D ( B) <+oo 


(25) 


(we denote by E the expectation with respect to the probability measure P). 
On another hand, since {Mjv(t)}t is a martingale, that is 

E(M N {t ) - -Mjv(s) /J’f) =0, VO <s<t<T 

and using the fact that 


Z t 

M N (t) - M n (s) = u N (t) - u N (s) - A{j u N (s)ds) + J B [ u N (s)g G (ujv(s))] ds 


has the same law as 


l Z 

M N (t) - Mn{s) =UN{t) -un{s) - A(j UN{s)ds)+ J B[d N (s)g G (ujv(s))] ds, 


it holds that 


E(M N (t) - M n {s) /F N s ) = 0, VO <s<t<T 
which implies that {Mpf(t)} t is a square integrable martingale on (A, J P) with 

_jV _ 

respect to T t = <j{un(s), s < t}. 

On another hand, since Mjv(t) has the quadratic variation (l24| . this means that 
for x, y £ D(B ), 


(M N (t),x)(M N (t),y)- I (a N (u N (z)),x)(a N (u N (z)),y)dz is an F? - martingale 
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Then 


E((M N (t),x)(M N (t), y) - (M N (s),x)(M N (s),y) 

t S 

- J(a N (u N (z)),x)(a N (u N (z)),y)dz + J{a N (u N (z)),x)(a N (u N (z)),y)dz /J 7 ™) = 0. 

0 0 

It follows that 

E((M N (t),x)(M N (t),y) - (M N (s),x)(M N (s),y) 

t S 

- / (aN(uN(z)),x)(aN(uN(z)),y)dz + / (aiv(wAr(2:)),x)(aiv(wAr(2:)),i/)(i2:/J’f) = 0. 


-=zN 


Henc e (MN(t),x)(M]sf(t),y)— / (o,n{un(z)), x)(cin{un(z)), y)dz is an T s —martingale 


and as a consequence, Mjy(t) has the unique quadratic variation process 


(M N (t),M N (t)) = J aN(uN{s))a* N (uN(s))ds. 
o 


From (l25l) . since 


ViV G N, sup E 

0 <t<T 


M N {t) 


D(B) 


sup E\M N (t)\ 2 D(B) < +oo, 

0 <t<T 


this means that (Mat)jvsn is a sequence of uniformly integrable martingales. 
Therefore, there exists on (A, T , P) a square integrable martingale M such that 
for N —> +oo, Mjsr —> M and 


t t 

M(t) = u(t) — uq — A(J u(s)ds) + J B [u(s)gG ( u ( s ))] ds 


with 


t t 

(M(t),M(t)}= lim / aN(uN(s))a* N (uN(s))ds = / A u + (s)ds. 
N^+ooJ J 

0 0 

By the representation theorem for square integrable martingales [S], there exists 
a cylindrical Wiener process W defined on (A, J 7 , Pte\o,T}i P) such that 


M(t) 


t 

J y/\u+{s)dW(s ), t G [0, T]. 
0 
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Then u satisfies that 


z z z 

i(t) = uo+A(J u(s)ds)—J B[u(s)go («(s))] ds+j \J Aii + (s)dl / l / (s), 


or equivalently 


u(t) = T(t)u 0 - [ T(t- s)B[u{s)g G {u{s))]ds 

t J 0 

+ f T(t- s)C(u{s))dW(s). 

J o 


Appendix A 

Lemma 7 For g £ [0, |], there exists constants C', C[ > 0, such that 

S 

J | T(t — u) — T(s — it)|^( B ) du < C' \t — 8\" 


and 

t 

J \ T (t - u)\ 2 d{b) du <C[\t- s\ 

S 

with 0 < s < t < T. 

Proof. We represen 

—w? to which correspond the eigenvectors • : 


Proof. We represent T(t) in terms of the eigenvalues of A = D-—r 

dx A 


T (t ) = ex P {~ w jt} (•> <Pj) <Pj 

j =i 


t€[0,T\ 


(26) 


(27) 

denoted 
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Then, we have 


S 

J \T{t~ u) -T(s-u)\ 2 D{B) du 
0 

S p oo 

< I |exp{—— «)} — exp {—u;|(s — w)}|^ (B) |<-,^)|^ (B) \^\ 2 D{B) du 
o j= : i 

S /» oo 

< / ^ |exp{-« 7 ?(s- w )}| 2 |l — exp{—Wj(t — s)}| 2 \(;Fj)\l {B) \<Pj\\ B) du 
0 J = 1 

oo S p 

^ D(B) \Vj\d(b) hf(i-s)| 2A / exp{-2u;?(a - u)}du 

i =1 n 


< sup |(.,^-)| D(B) sup |^. I D(B) 1 1 - s) | 2A 

J j=l 


ID(B)' 

OO 


2A —(2—4A) 

x V 


< Cb|t-«)| 2A X) 

i=i 


2A —(2—4A) 

' V 


where we have used the fact that for all u € [0,1], 

1 — exp{— x} <x u (x > 0 ). 

OO 

For A € [0, j], we have Wj *■ 2 4A ' ) < oo. If we set rj = 2A fq £ [0, |]), hence 

3= 1 
s 

J \ T (t -u)- T(s - u)\ 2 d(b) du <C'\t — s)| I? 

0 

Similarly, for the same rj £ [0, as above, we have 

t ^ t 

A OO A 

J J exp{- 2 «;?(t - u )}du 

s J = 1 a 

1 X 

< -sup |(.,^-)| 2 sup \Tj\ 2 d(b) ^2w~ 2 (1 - exp {-2 w](t - a)}) 


3 =1 


< c'o —«r 


w. 


-(2—277) 


<Cl|i-a)r 


i=i 


Lemma 8 For 77 € [0, 5 ], 1 /iere exists constants C",C 1 > 0 ,such. that 
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( 28 ) 


J |B [T(t -u)- T{s - «)]\ 2 d(b) du < C" 1 1 - s'" 1 


Z 

J | BT{t - u)\ 2 d(b) du < C'[ 1 1 - s| 


with 0 < s 


Proof. From the representation 


we obtain 


Hence 


Tit) = Yl ex P{- w ‘jt} <Pj- 

3 =1 


BT(t) = '*rex.p{-w?t}(.,<p j )B<p :j . 


S 

J I B [T(t -u)- T(s - u)]\ 2 d{b) du 


-w](t - u)} - exp{—ru|(s - «)}|^ (B) \B<P^ d{b) du 


-w 2 {s - u)}| 2 |1 - exp - s )}| 2 \ 2 d(b) | B<pjf D(B) du 


OO s » 

^ S D(B) \ B Tj\ 2 D{B ) I w j(t-s)\ 2X J exp{—2w|(s — u)}dh 


< sup |(., <pj)\ 2 sup \B^\ 2 1 1 - s)| 2A 


2A —(2—4A) 

2^ W 3 


<C 2 \t-s)\ 2X Yl 

3=1 


2A —(2—4A) 

2^ w j 


where we have used the fact that for all u £ [0,1], 

1 — exp{—a;} < x u [x > 0). 

OO 

For A £ [0, |], we have w -( 2 - 4 A) < oo. If we set 77 = 2A (77 £ [0, |]), hence 


S 

J I B [T(t -u)- T(s - u)}\ 2 d{b) du < Ci \t - s)\' 
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Similarly, for the same rj £ [0, as above, we have 


Z „ Z 

/ OO p 

\ B T(t - u)\l (B )du <J2\(-’Vj)\ 2 D (B)\ B Vj\D(B) / exp {-2 w?{t-u)}dv 

7 = 1 J 

s J s 

OO 

< isup |(.,^)|" sup ^ W t 2 (1 - exp{- 2 W 2 (f - a)}) 


i=i 


<C'\t-s\^ 


(2—277) 


<c;it- S )i 


j=i 

V 


Lemma 9 (Garsia-Rodemich-Rumsey) Let the function if : [0, oo[—» [0, oo[ 
non-decreasing with lim if(u) = +oo. Let the function p : [0,1] —> [0,1] be 

u —>+oo 

continuous and non-decreasing with p( 0 ) = 0 . 

Set < ip^iu) = supu if if if)) < u < oo . Let f be a continuous function on [0,1] 

[ tp(v)<.u 

and suppose that 


^ l/(aQ-/(y)l 

p(x - y) 

o o 

Then, for all x,y € [0,1], we have 



)dxdy < B < oo 


\y-x\ 

I/O) - f(y)\ < 8 J if~ l (^)dp(u). 

o 


Appendix B 
Proof B 1 Since 


I fit) — fit') | D , m 

ll/llc*([0,T],D(B)) = S ‘ UP D(B) + Su P - ---f-' 

te[o,T] \t-t'f 

this means that V t,t' £ [0, T],t ^ t' 

I/O ~ /(^ )l d(b) — II/IIc t ([o,t],d(b)) l*-*T 

As f £ Br, we have 


\m-m\D {B) <R\t-t'\ s 

This means that Br is uniformly equicontinuous. By the theorem of Ascoli- 
Arzela, Br is a compact of C([0,T],D(B)). 
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